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Abstract 

A requirement for many quantum computation schemes is the abil- 
ity to measure single spins. This paper examines one proposed scheme: 
magnetic resonance force microscopy, including the effects of thermal 
noise and back-action from monitoring. We derive a simplified equa- 
tion using the adiabatic approximation, and produce a stochastic pure 
state unraveling which is useful for numerical simulations. 



1 Introduction 

Single-spin measurement is an extremely important challenge, and necessary 
for the future successful development of several recent spin-based proposals 
for quantum information processing. [1, 2, 3, 4, 5] There are both direct and 



indirect single-spin measurement proposals. The idea behind some indirect 
proposals is to transform the problem of detecting a single spin into the 
task of measuring charge transport [2, 6], since the abihty to detect a single 
charge is now available. For direct single-spin detection, magnetic resonance 
force microscopy (MRFM) has been suggested [7, 8, 9] as one of the most 
promising techniques. To date, the MRFM technique has been demonstrated 
with sensitivity to a few hundred spins [10, 11]. 

In this paper we discuss how to read out the quantum state of a single 
spin using the MRFM technique based on cyclic adiabatic inversion (CAI). 
[12, 10, 9] In this CAI MRFM technique, the frequency of the spin inversion 
in the rotating frame is in resonance with the mechanical vibration of an 
ultra thin cantilever, allowing it to amplify the otherwise extremely weak 
force due to the spin. These amphfied vibrations can then be detected by, 
e.g., optical methods. 

Previous studies [8, 9] of the dynamics of single-spin measurement by 
MRFM considered only the unitary evolution of the spin and the cantilever 
system, without including any effects of external environments or measure- 
ment devices. Only recently, the effect of thermal noise environment on the 
dynamics of the spin-cantilever system in the MRFM was studied [13] by us- 
ing the Caldeira-Leggett master equation [14] in the high-temperature limit. 

There is, however, a macroscopic device in the MRFM setup which mea- 
sures the cantilever motion and hence provides information about the spin 
state. To our knowledge, the back-action of the measurement device and the 
effect of the thermal noise on the dynamics of the cantilever-spin system for 
the single-spin detection problem by MRFM have not yet been investigated 
systematically. In this paper, we include, in our analysis, a measurement de- 
vice (a fiber-optic interferometer) to monitor the position of the cantilever. 
We consider various relevant sources of noise and calculate the signal-to-noise 
ratio of the output photocurrent of the measurement device. We also develop 
a realistic continuous measurement model, and discuss the approximations 
and conditions to achieve a quantum non-demolition measurement of a single 
spin by MRFM. Finally, we present some simulation results of the dynamics 
of the single-spin measurement process. 
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Figure 1: Schematic diagram of the MRFM setup. 



2 The measurement scheme 

A schematic illustration of the MRFM setup is shown in Fig 1. A uniform 
magnetic field, Bq, points in the positive ^-direction. A single spin is placed 
in front of the cantilever tip which can oscillate only in the 2;-direction. A 
ferromagnetic particle (or small magnetic material) mounted on the can- 
tilever tip produces a non-uniform magnetic field or magnetic field gradient 
of (dBjdZ) on the single spin. As a result, a reactive force (or interaction) 
acts back on the magnetic cantilever tip in the z-direction from the single 
spin. The origin is chosen to be the equilibrium position of the cantilever tip 
without the presence of the spin. 

In CAI, the cantilever is driven at its resonance frequency to amplify the 
otherwise very small vibrational amplitude. This is achieved by a modulation 
scheme using the frequency modulation of a rotating radio-frequency (RF) 
magnetic field in the x-y plane. In this case, the rotating RF field can be 
represented as Bi^ = Bicos{[uj + Auj(t)]t}, Biy = —Bism{[u + Au{t)]t}, 
where the frequency modulation Aooit) is a periodic function in time with 
the resonant frequency Um of the cantilever. In the reference frame rotating 
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with the Bi, the spin-cantilever Hamiltonian can be written as 



(dB \ 
~dFj ' ^^"^ 

where ujl = dfJ^Bz/h and ui = giiBi/Ti are respectively the Larmor and Rabi 
frequencies, B^ includes the uniform magnetic field Bq and the magnetic 
field produced by the ferromagnetic particle, g and // are the g^-factor and 
the electron or nuclear magneton, and 

is the Hamiltonian of the cantilever in isolation (i.e., with no external mag- 
netic field coupling it to the spin). For uj — ujl, we arrive at an effective 
cantilever-spin Hamiltonian of the form 

Hsz{t) =Hz- 2r]ZS, + f{t)S, - eS^ , (3) 

where f{t) = Auj{t), rj = {g^/2){dBz/ dZ)^ and e = hui. We will discuss 
in details the rotating picture and adiabatic approximation for the spin- 
cantilever system in the next section. 

In the following, we briefiy describe the basic principle of the single-spin 
measurement by CAI MRFM. In the case when the adiabatic approximation 
is exact, the instantaneous eigenstates of the spin Hamiltonian in the rotating 
frame of the Bi field are the spin states parallel or antiparallcl to the direction 
of the effective magnetic field B'^^(t) = {£,0, —f{t)), denoted as \v±{t)), 
respectively. We define an operator S'^ for the component of spin along this 
axis. Note that the initial spin state in the laboratory frame has the same 
expression as the initial state in the rotating frame. Starting at a general 
initial spin state (in the laboratory or rotating frame) of 

X(0) = a| T) + &U) (4) 

in the Sz representation, we can rewrite this initial state in the basis of the 
instantaneous eigenstates of S'^ as 

x(0) = aeff|^;+(0))+6eff|^-(0)), (5) 
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where 



Oeflf = acos(©o/2) + &sm(©o/2), 
= -asin(eo/2) + 6cos(©o/2), 



(6) 
(7) 



and 00 = 0(0) is the initial angle between the effective magnetic magnetic 
field and the z-axis direction. This implies tan[0(t)] = Bf{t)/Bf{t) = 
—e/f{t). It then follows from the adiabatic theorem that the spin state at 
time t can be written as: 



where X±{t) are instantaneous eigenvalues. So the probabilities of finding the 
spin to be in the instantaneous eigenstates \v±{t)) are respectively la^sl^ and 
l^eflfP- Since the coefficients and b^s are time independent, the probabil- 
ities loeffp and |6efrP remain the same at all times. This provides us with an 
opportunity to measure the initial spin state probabilities at later times. 

How do we measure these spin state probabilities? The idea is to transfer 
the information of the spin state to the state of the driven cantilever. In the 
interaction picture in which the state is rotating with the instantaneous eigen- 
states of the spin Hamiltonian, the spin-cantilever interaction can be written 
as 2rjZS'^cos[Q{t)]. As a result, the phase of the driven cantilever vibrations 
depends on the orientation of the spin states. Suppose that the initial state 
is a product state of the cantilever and spin parts. At a later time, due to 
the interaction between them, the total state becomes entangled. Monitoring 
the phase of the cantilever vibrations will give us the information about the 
spin. Numerical simulations (see Fig. 4 with reasonable parameters for the 
CAl approximations) indicate that as the amplitude of the cantilever vibra- 
tions increases with time, the phase difference in the oscillations for the two 
different initial spin eigenstates of S'^ approaches tt. In other words, the mea- 
surement of the single-spin states can be achieved by monitoring the phases 
of the cantilever vibrations at some later time t. Phase-sensitive, optical ho- 
modyne measurements of the cantilever vibrations can be performed using 
a fiber-optic interferometer. The main purpose of this paper is to present 
a realistic and detailed analysis of the single-spin measurement scheme, in- 
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eluding the effects of the measurement device and other relevant sources of 
noise. 



3 The rotating picture and adiabatic approx- 
imation 

We assume an effective cantilever-spin Hamiltonian of the form (3) where for 
the moment we let f{t) and e be arbitrary, and Hz is the Hamiltonian given 
by (2). It is useful to group this into three terms 

Hsz{t)^Hz + Hi + Hsit) , (9) 

where 

Hj = -2riZS,, 
Hs{t) = f{t)S,-eS,. (10) 

The state of the cantilever-spin system evolves according to the Schrodinger 
equation 

^^4h..WI.W>. (u) 

In realistic cases, the spin part of the Hamiltonian (representing preces- 
sion under the magnetic field) gives an evolution which is very rapid compared 
to the reaction time of the cantilever. It therefore makes sense to switch to an 
interaction picture in which the state is rotating along with this precession. 
We do this by introducing a (partial) time translation operator 



Us{t) =: exp-- 







Hs{t')dt' 



(12) 



where :: indicates that the integral is to be taken in a time-ordered sense; 
this unitary operator obeys the differential equation 

We then introduce the state \il^{t)) in the rotating picture: 

m))^ulm{t)), (14) 
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with \ip{t)) the solution of the original Schrodinger equation (11) at time t. 
The evolution equation for \ip{t)) is 

= ^uimsmit)) - ^uimszmit)) 

= -^Hzm)) - { [uimMt)] \m) 

= ~Hz\m) + [iJi{t)s,iJs{t)\ m)) ■ (15) 

We can define a locked spin operator 5"^ {t) 

SL{t) = [ul{t)S,Us{t)] ; (16) 
in terms of this, the equation of motion for 1-0) becomes 

'-^ = -^H,\m + 'i^zU)\m)- (17) 

Unfortunately, it is difficult to get an exact solution for Usit) for a general 
function f{t). This means that it is also difficult to derive an exact expression 
for Slit), and the rotating picture (15), while formally correct, is not very 
helpful. 

However, while we cannot easily find an exact expression for tJs{t) for 
general f{t), we can easily find an approximate solution for a large class 
of functions. Suppose that e is large and f{t) is slowly varying, so that 
|/(t)|,e > \f'it)/f{t)\ for typical values of /(t) and f'(t). Then Hsit) is also 
slowly varying, and if a spin begins in an instantaneous eigenstate of Hs{t), 
it will remain close to an instantaneous eigenstate of Hs{t) for all times by 
the adiabatic theorem. 

The instantaneous eigenstates of Hsit) are 

Hs(t)\v±(t)) = X±(t)\v±(t)) = ±X(t)\v±(t)) , (18) 

where 



-2 



Kit)) = , ' ID- , ^^'^^'^'^ IT). (19) 

/{f{t)T\{t))' + e' J{f{t)TX{t))' + e' 
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We use these instantaneous eigenvectors and eigenvalues to define an approx- 
imation to the unitary operator Us{t): 

tj'sit) = / ® |^;+(i))K(0)|e-**W + / |^;_(i))(^;_(0)|e**W , (20) 

with the accumulated phase 

m^lJ^mdt'. (21) 

Note that $(t) obeys d^{t)/dt = X{t). This implies that 

= -^A(t)/®|t;4t))(t;40)|e-^*W + lA(t)/®|t;_(t))(t;_(0)|e^*W 
+i ^^^^K(0)|e-^*W + / %^(^;_(0)|e^*W 

- -^^5(t)?7^(t) + /^^^K(0)|e-*« 

+/®^^^^(^-(0)|e^*^*), (22) 

which has the form of (13) plus some additional terms. From the definition 
(19) of \v±{t)), we see 

Provided that f{t) is slowly varying, the additional terms in (22) will be 
small. 

Just as before, we can define a rotating picture, now using the unitary 
transformation U'sit), 

\m ^ {u'sif))' \m) ■ (24) 

This gives us a new evolution equation for 1-0): 

dt dt dt 



+l--(0))^^e^*w)t/^(t)|^(t)). (25) 
At this point, it is helpful to introduce a new set of spin operators 

s'y = |/®(K(o))Mo)|-Mo))(^-(o)|) , 

s'. = ^/®(K(o))K(o)|-K(o))(^_(o)|) . (26) 

Using the definition (20) for Us{t), we can solve for the various terms in (25): 
{Um'S.U'sit) = -^Si - -|- (^icos(2$(t)) - ^;sin(2$(t))) . (27) 



k^(0))%^e-*(*) + |._(0))%^e-« ) U'sit) 



ie df{t) 



(^;sin(2$(i)) + ^;cos(2$(i))) . (28) 



nx^it) dt 

Substituting (26-28) into (25), we get 

+T^^A^ (^ico«(2^W) - '5;sin(2$(t))) m) 

^'Wif)^ (5;sin(2$(t)) +5;cos(2$(t))) |^(t)) . (29) 

Note that this equation is still exact — it is equivalent to the original 
Schrodinger equation (11). However, we can see that if \ f{t)\,e are large, then 
$(t) will be a rapidly growing function, and the last two terms of equation 
(29) will oscillate very rapidly compared to the first two terms. Over a short 
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period relative to the response time of the cantilever they will essentially 
average away to nothing. In this limit, therefore, we can reasonably make a 
rotating-wave approximation, to get the approximate evolution equation 

^ {Hz - 2r^{f{t)/mZS',) m)) . (30) 

This is equivalent to making an exact adiabatic approximation, as described 
in section 2. We can see how this approximation compares to the complete 
Hamiltonian for a reasonable set of parameter values in figures 2 and 3. This 
set of parameters was chosen to match those of Berman et al. [9] , as was the 
set of times plotted in figure 3. Comparison shows that our results match 
their unitary simulations to good precision. For the rest of this paper we 
will be using the rotating wave approximation, and representing states in 
the rotating frame. For simplicity, we henceforth omit the accent from the 
state 




Figure 2: Mean cantilever position (Z) vs. t for the complete and rotating 
wave Hamiltonians. 

In this rotating-wave approximation, if the spin begins in an instanta- 
neous eigenstate of Hs{t), it will remain in an instantaneous eigenstate at all 
times. If it begins in a superposition of the two eigenstates, the spin and can- 
tilever degrees of freedom will become entangled, with the two components of 
the wavefunction corresponding to the two spin directions remaining undis- 
turbed for all times. Monitoring the position of the cantilever then serves as 
a nondemolition measurement of the spin, as we would wish. 
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Figure 3: The probability distribution, of finding the cantilever at posi- 
tion 2; at a range of times for the complete and rotating wave Hamiltonians. 



11 



Note that the corrections to the adiabatic approximation include terms 
which can flip the spin. These terms must remain small for the system to 
be a true nondemolition measurement. The result of the spin measurement 
manifests itself as a vr phase shift in the oscillation of the cantilever. We can 
see this in figure 4. 
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Figure 4: Mean cantilever position (Z) vs. t for initial spin up and down in 
the S'^ direction. 
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4 The thermal environment 



Unfortunately, in practice we cannot treat the cantilever as an isolated sys- 
tem. It is coupled at least weakly to the vibrational modes of the bulk, and 
is therefore subject to dissipation and thermal noise. Since the cantilever 
can be treated as a single harmonic oscillator, we can model the effects of 
this thermal bath by the well-known Caldeira-Leggett [14] master equation 
in the high-temperature limit: 

p-~[Hszit),p] - '^[Z,{p,p}] - ^,[Z,[Z,p]] , (31) 



where the parameters are 



r 

h 



e = ^ , (32) 

m is the cantilever mass, T is the temperature, k is Boltzmann's constant (or 
the equivalent for our system of units), and F is the strength of the coupling 
to the thermal bath. We can interpret 7^ (with units of inverse time) as 
the dissipation rate and i (with units of length) as the thermal de Broglie 
wavelength. 

A feature of this equation is that it doesn't necessarily preserve the posi- 
tivity of p on short time scales (though at long times it is well-behaved) [17]. 
This arises because of the approximations which are made in the derivation, 
which become invalid at very short times. While this may be physically 
unimportant, it can be inconvenient; in particular, if we wish to unravel the 
evolution into a stochastic Schrodinger equation [23] (as we will in section 
6), it is necessary to start with a master equation in Lindblad form, [16] 

p = ~[H, p]+Y: [2L,pL] - {LjL,, p}] , (33) 
^ j 

for some Hermitian H and set of general Lindblad operators {Lj}. The 
Caldeira-Leggett equation (31) is not of this form, which is why it can violate 
positivity of p. 

The exact quantum Brownian motion master equation was shown [17] not 
to have the Lindblad form, but rather requires time-dependent coefficients 
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to ensure the positivity of the density matrix at short times. However, by 
keeping more terms from the high- or medium-temperature-hmit expansion 
in a consistent way, Diosi [18] has shown that the Caldeira-Leggett equation 
can be replaced by another master equation which is of Lindblad form, and 
which agrees with it except at very short times when the equation's vahdity 
is questionable in any case. This is done by adding a term to (31) of the form 
— (7to£^/2/i^)[p, [p, p]]. The procedure is analogous to completing the square. 
If we choose the ansatz 

L^AZ + iBp (34) 

with real A, B, plug it into the equation (33), and equate it to the Caldeira- 
Leggett equation (31) plus the additional term, we get 

p = ~m)[H,p]~A'[Z,[Z,p]]~B%[p,p]]^ 

+iAB{—2Zpp + Zpp + pZp + 2ppZ — pZp — ppZ) 

= -{i/h)[Hsz{t),p] - ^[Z, [Z,p]] - \p,p]] 

+ ^ {PpZ - Zpp - Zpp + ppZ) , (35) 

which implies that 

A = ^-fj2P , 
B = yfiJ^\ 

H = Hsz{t) + i^j2){Zp + pZ) = H'szit). (36) 
So the Lindblad operator for this equation is 

L = ,fij2 {{i/e)z + i{e/h)p) , (37) 

and the effective Hamiltonian, going to the rotating picture and making use 
of the approximation derived in section 3, is 

^'''^^^ ^ 2^" + " 2r^{f{t)/\{t))ZS', + {^m/2){Zp-rpZ) . (38) 

In order for the cantilever to be an effective measurement device, the loss 
rate must be very low: a;^ ^ j^- 
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5 The effects of monitoring 



In order to serve as a measurement scheme, we must have some way of 
monitoring the motion of the cantilever. Because of the microscopic scale 
of the motion, this is not so easily done. One approach is to use optical 
interferometry to measure the cantilever position. 

As shown in figure 1, the cantilever forms one side of an optical microcav- 
ity and the cleaved end of the fiber forms the other side. As the cantilever 
moves, the resonant frequency of the cavity changes. Because the timescale 
of the cantilever's motion is very long compared to the optical timescale, we 
can treat the effects of this in the adiabatic limit. The cavity mode is also 
subject to driving by an external laser, and has a very high loss rate. The full 
master equation [19] for the cantilever-spin-cavity system in the interaction 
picture is 



where H'gzit) and L are the Hamiltonian and Lindblad operator for the 
cantilever and spin given by eqns. (37) and (38), E is the strength of the 
laser driving. A is the detuning from the "neutral" cavity frequency, k is the 
coupling strength of the cantilever to the cavity mode, and 7c is the loss rate 
of the cavity. 

Suppose now that we perform homodyne measurement [15, 20] on the 
light which escapes from the cavity. We would like to replace the equation 
(39) above with an equation for the conditional evolution of p, conditioned 
on the output photocurrent Ic{t)- The conditional evolution equation for our 
system then becomes [20, 21] (in Ito calculus form) 



P 



-^[Hszit),p] + 2LpD - Dip - pDL 

-i[E{a} + a) + a^a(A + kZ),p\ 
-|-(7c/2)(2apa^ — a' dp — pa)d) , 



(39) 



dp = ~[H'sz{t),p]dt + {2LpD - Dip - pDL) dt 
-i[E{d^ + a) + a^a(A + kZ), p]dt 




(40) 
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where < e^i < 1 is the detector efficiency and dWt is a real stochastic 
differential variable which obeys the statistics 



M[dWt] = , M[dWtdWs] = d{t - s)dt . 



(41) 



This noise is related to the output photocurrent [15, 20, 21] 



dWt 



(42) 



dt 



where (3 is a constant giving the device's range of response. 

We want to operate in the "bad cavity" limit where jc 3> uJm- This 
means that the cavity mode will approach equilibrium on a timescale very 
short compared to that of the cantilever's motion, so that the cavity mode 
can be adiabatically eliminated [19, 20, 21] from this equation, leaving an 
equation in terms of the spin and cantilever position alone. 

Let the detuning vanish A — > and the coupling k. to the cantilever be 
very small. If we initially neglect this coupling altogether, we can solve for 
the steady-state of the cavity mode in isolation from the cantilever: 



which implies that p — \ao){cKo\, where d\ao) — CKo\ao) is a coherent state 
with 



Now let us restore the coupling k between the cantilever and the cavity 
mode. If this coupling is very small, then the state of the cavity mode will 
remain very close to the state |q;o). In this case, it is very useful to switch 
to a displaced basis [19, 20, 21] for the cavity mode. We switch from the 
operators a, d'^ to displaced operators 




(43) 



b = d — ao , 



6^ = — ctQ , 



(45) 



and displaced number states 




(46) 
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Obviously |0) = |«o), and |1) = a^lcto) — ctglao)- 

We now make the ansatz of keeping the two lowest displaced number 
states |0, 1) of the cavity mode and neglecting the rest. [20, 19, 21] We then 
write the full density matrix for the spin-cantilever-cavity system as 

p{t) = Po{t) ® |0)(0| + ® |1)(0| + p\{t) ® |0)(1| + p2{t) ® , (47) 

where Po,i,2 ^-re operators which act on the Hilbcrt space of the cantilever 
and spin, and /7o,2 are self-adjoint. The reduced density matrix of the spin- 
cantilever alone is obtained by tracing out the cavity mode, yielding 

Psz{t) ^ Po{t) + p2{t) . (48) 

If we substitute the definitions (45) and (47) into the stochastic master equa- 
tion (40) and collect terms, we get a set of coupled equations in the operators 

dpo = (^-^[H'sz{t),Po]+'^LpoL^ -L^Lpo-poL'^L^dt 
■[Z,po\dt H (Zpi -I- p[Z)dt + 7cP2«^ 



7c " ■ ' 7c 



+V^d (pi + p\ - PoTr{pi + p\ }) dWt , (49) 

dpi = (^-^[H'sz{t),Pi] + 2Lp^P - L^Lp^- p^PL^dt 

—inZpidt — [Z,pi]dt {Zpo — P2Z)dt — (7c/2)pi(ii 

7c 7c 

+V7c^ (P2 - PiTr{pi + p\}) dWt , (50) 



dp2 = [-^[H'sz{t),p2]+2Lp2D -PLp2-p2L^Ljdt 

( 4:iKE^\ . , 2kE.^ + , 

- \iK-\ — [Z, p2\dt {Zp[ + piZjdt - ^cP2dt 

V 7c / 7c 

-V7^P2Tr{pi + pl}dWt . (51) 
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Both pi and p2 contain damping terms, which imply that they will remain 
small at all times provided nZ is sufficiently small compared to 7c. (This also 
implies that our ansatz is reasonable for sufficiently small k.) 

By making use of the above equations, we can find the evolution equation 
for the reduced density matrix psz'- 

dpsz{t) = dpo{t) + dp2{t) 

= [-\[H'sz{t).Psz] + 2LpszL^ - L^Lpsz - PszL^L^ dt 

■ [Z, psz\dt H [Z, pi - p[\dt - ik[Z, p2\dt 



+V^d (pi + p\ - PszTr{p, + p\}) dWt . (52) 

If we keep only terms to second order in K.Z we can neglect the p2 term. 
This leaves only the terms proportional to pi ± p| , which we need only know 
to leading order in kZ. Provided (as we have already assumed) that the 
cantilever moves slowly compared to the timescale set by 7^ and that kZ 
can be treated as small, then to leading order dpi vanishes; pi remains in 
an approximate equilibrium state. If we make use of this assumption we can 
(again to leading order) solve for pi± p\: 

+ 4:KE , ^ -, 

Pi + Pi ~ ^{Z,psz}: 

Tc 

Pi- Pi ~ ^[Z,psz\, (53) 

7c 

which when inserted into (52) gives us a closed evolution equation for psz'- 

dpszit) = [-^[H'szit),Psz] + ^LpszL^ - t'^Lpsz - PszL^L^ dt 
■[Z, psz\dt 5— [Z, Psz\\dt 



il it 



+VT^d^ {Zpsz + PszZ - 2p5zTV{p52}) dWt . (54) 

7c 

(Note that we have absorbed a factor of —1 into dWt-) 
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Examining the terms in (54), we see that by ehminating the cavity mode 
we get another effective term in the Hamiltonian, and another Lindblad 
operator. We can therefore write this stochastic master equation in the form 

2 

dpsz{t) = -kHe{{it),Psz]dt + J2 [^LjPszL] - L]LjPsz - PszL]Lj) dt 

((L2 - {L2))psz + Psz{L2 - (L2))) dWt , (55) 
where we now make the definitions 

+^-^+ (7m/2)(^p + p^) . (56) 
7c 

Note that the term AnE'^Z/'f'^ is a constant force, which just displaces the 
equihbrium position of the cantilever. It can be eliminated simply by chang- 
ing the origin of Z, and is in any case small for reasonable values of the 
parameters. The output from the homodyne measurement now corresponds 
to a measurement of the cantilever position (Z): 

W^f^l — \^/ + ^^^"^)- v57) 

As we shall see in the next section, we can further unravel this stochastic 
master equation (55) into a stochastic Schrddinger equation for pure states. 
This further unraveling provides considerable improvement in numerical ef- 
ficiency, though it does not represent an actual measurement process. 

6 Pure state unraveling 

The stochastic master equation (55) represents the evolution of the cantilever- 
spin system, conditioned on the photocurrent measurement record Ic{t). If 
we averaged over all possible measurement records, the dWt terms would 
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average to zero, and we would be left with an ordinary deterministic master 
equation for the cantilever and spin. It is for this reason that the stochastic 
master equation is therefore often referred to as an unraveling of the average 
master equation. 

For numerical purposes, it is often much easier to solve an equation for a 
pure state vector raihci than a density matrix [22, 23]. It is therefore useful 
to unravel equation (55) still further to an equation which preserves pure 
states. We do this by introducing a second stochastic process. 

First, let us ideahze to perfect detector efficiency Sd — 1. We then intro- 
duce the new master equation 

dpsz{t) = -T[H<,^{t),psz\dt + hLjpszL] - L]Ljpsz - pszL]Lj\ dt 

+ ^ ((4- - {L,))psz + Psz{L, - (4))) dW,t , (58) 

where the Hamiltonian and Lindblad operators are the same as in (56) and 
we now have two independent noise processes represented by stochastic dif- 
ferential variables dWu and dW2t which satisfy 

M[dWjt] = , dWitdWjs = S{t - s)6ijdt . (59) 

If we take the mean of (58) over dW^ we recover equation (55). We can think 
of the additional stochastic process as representing a fictitious additional 
measurement, whose outcome we average over to recover the state which is 
conditioned on the actual measurement. 

However, equation (58) has a great advantage over (55). If psz is initially 
a pure state psz = \ipsz){'4'sz\, it will remain a pure state at all times, the 
state of course depending on the stochastic processes Wi and W2- We can 
recover the solution of (55) by averaging 

psz{t) = MwMsz{t)){^szm . (60) 

It would be useful to replace equation (58) with an explicit evolution equa- 
tion for \ipsz) instead of psz- This equation is the quantum state diffusion 
equation with real noise [24, 25]: 

dli^Sz) = -^^eff(i)|V'5z)rft + ^E (2(4^)4 1(4)1') IV'5z)rfi 



20 



+ ^T.{Lj-{Lj))\i^sz)dW,t 



(61) 



The nonlinearity of this equation arises to preserve the norm. 

7 Numerical simulation 

We have simulated this system using the C++ quantum state diffusion h- 
brary [26] to numerically solve both the unitary evolution with Hamiltonian 
(30) and the stochastic equation (61). All of the figures in this paper were 
generated using this software. 

We chose our parameters based on those used by Berman et al. , [9] . These 
values are (in arbitrary units): 



h 




V 


= 0.3, 


e 


= 400.0, 




= ^m/Q = 0.00001 , 




= 10000.0, 



(62) 



where Q is the quality factor of the cantilever. The driving force f{t) takes 
the form 

. . _ r -6000 + aoot if o < t < 20; , . 

•'^^^ ~ U000sin(t-20) ift>20. ^ > 

If we make contact with physical values for actual cantilevers used in experi- 
ments, we have o;^ fa 10^s~^ and m ^ lO^^^kg. The value of ksT above then 
corresponds to a temperature of around O.IK, which is within the bounds 
of experimental feasibility, though rather lower than the temperatures used 
in current experiments (around 3K) [11]. Since rj = {gfi/2){dBz/dZ)Q, the 
value of rj corresponds to a field gradient of about 1.5 x lO^T/m, which is 
higher than current experiments by roughly two orders of magnitude [11], 
but hopefully this too will improve with time. The cantilever would undergo 
displacements of about a nanometer. 

Alternatively, rather than increasing the field gradient we could achieve 
similar numbers by lowering the spring constant of the cantilever, for instance 
by shrinking the mass of the cantilever. Lowering the mass by a factor of 100 
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has the same relative effect on rj as increasing the field gradient by a factor 
of ten. 

We then might ask about realistic parameters for the monitoring. A 

typical cavity size L is about a micrometer, with a laser frequency of oOc ~ 
1.4 X IQ-'^^s^^. This cavity is generally quite lossy; reasonable quality factors 
might be in the range Qc ~10-100. The parameter E is a function of the 
laser power, E = ^ Pjc/ ^^^c = \J P/^Qc- For P ~ 1//W and Qc ~ 100 we 
have E ~ lO^^s"^. The coupling between the cantilever and the cavity is 
given by a geometric factor k — Uc/L ~ 1.4 x 10^^(m • s)~^. In arbitrary 
units, this gives coefficients 



%kE 

7c 
AkE"^ 



1.9 X 10 



7 X 10^ 



3 



\ 7 



= 0.07. (64) 



The first value is the multiplier in (57); the second gives the equilibrium 
displacement of the cantilever; the third is the coefficient of the Lindblad 
operator L2. 

One question we can now easily address is how quickly the state of the 
spin collapses onto eigenstates of S'^. In figure 5 we plot {S'^ for ten different 
trajectories. We see that in all ten cases the spin converged to ±1/2 quite 
quickly, before t = 80. 

If we compare this to the results of figure 4, we see that the spin state 
collapses rather more quickly than the cantilever oscillations can respond. We 
only get a clear output signal when the two phases are well separated, which 
does not occur until nearly t — 150. Generically, the difficulty of collapsing 
the spin state is much less than the difficulty of obtaining an unequivocal 
readout. 

The curves depicted in figure 4 are idealized, without the measurement 
noise which will always be present in the output current (42) or (57). In figure 
6 we show what actual output would look like for the set of parameters we 
are discussing. Note that even with the noise, the two phases (representing 
spin up and spin down) are clearly distinguishable. In the next section, we 
derive an expression for the signal-to-noise ratio in more general situations. 
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20 40 60 80 100 



Figure 5: Expectation value (S'^) vs. t in arbitrary units for ten different 
trajectories, sliowing the rapid localization of the spin, for an initial super- 
position state (|v+(0)) + |w-(0))) /V2. 



8 Signal-to-noise ratio 

Since we have to detect the effect of very weak force on the cantilever by the 
single spin, we need very high resolution for the cantilever position measure- 
ments and a good control of the various noise sources in the MRFM device. 
As described in section 2, the tiny displacement of the cantilever is measured 
by a fiber-optic interferometer as a phase shift of the interference fringes. We 
shall analyze the quantum and thermal noise in this homodyne measurement 
scheme. 

The Hamiltonian for the combined system of the spin, cantilever and 
cavity mode, excluding coupling to the environments, in the spin rotating 
frame is 

' X{t) ' 

+;iE(a^e-*'^°* + ae^'^«*) + fiKa^aZ. (65) 

Here Uc is the optical frequency of the cavity mode, Uq ~ Uc is the driving 
frequency of the external laser and other terms and parameters have been 
described in section 5. The master equation approach in section 4 is valid in 
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Figure 6: Simulation of photocurrent output in arbitrary units, including 
measurement noise, using parameters of section 7. We have chosen the scale 
(3 so that the vertical scale matches that of figure 4, and also plotted the 
signals without the noisy dW/ dt components. 
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high or medium temperature case. Here we analyze the noise in the Heisen- 
berg picture, using the quantum Langevin equation approach that is vahd at 
any temperature. [27] 

Using standard techniques, [28, 29] the reservoir (environmental) variables 
may be eliminated, in the interaction picture with respect to hujQa'^a, to 
give the quantum Langevin equations describing the dynamics of the whole 
system: 

at m 

^ = -ma;^Z(i)-Ip(i)-W(t)a(i) + W(i) + 2r7(^^^(i)(67) 
at m X\t) 

^"^^^ - -{iuJc-iujQ + ^)a{t) -mZ{t)d{t) -iE + ^^tti^^t)^ (68) 



dt ^ " " 2 

dm 



dt 



0, (69) 



In the equations, the usual optical input noise operator di^{t) is associated 
with the vacuum fluctuations of the continuum of electromagnetic modes 
outside the cavity and its correlation function is given by 

{h,,{t)a}^^{t'))=6{t-t'). (72) 

The random force VV(t) describes the thermal noise motion (quantum Brow- 
nian motion) of the cantilever at temperature T. For the case of an Ohmic 
environment, the thermal random force correlation is given by [27] 

{W(t)W{t')) = —[Tr{t - t') + iTi{t - 1% (73) 

TT 



where 



^rit) = / (ici;ci;cos(a;t)coth — , (74) 

Jo ylksT J 

^i(t) — / duj uj sm(ujt) , (75) 

^0 
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with the frequency cutoff of the reservoir spectrum. Without the presence 
of the external driving force from the spin, the cantilever-cavity system can be 
characterized by a semi-classical steady state with a new equilibrium position 
for the cantilever, displaced by Zgt = —ti\ast\^/{rnuj'^ with respect to that 
with no external driving laser field, and the cavity mode in a coherent state 
\ast) with the amplitude given by 

-^E 

where A = ujc — ojq — K^\ast\^ / {muj'^ is the cavity mode detuning. By 
adjusting either ujq or ojc, the detuning can be set to zero A = 0. As a result, 
Ust = «o = —'^iE/'-jc- Linearizing the quantum Langevin equations about 
the steady-state values and renaming with Z{t)^ a{t) the operators describing 
the quantum fiuctuations around the classical steady state, we obtain 

^ - >■ 

dpit) - r 

— - — = —muj'LZit) pit) — hK[aQaHt) -\- aQait)] 

dt m 

+W(t)+2r;{||-S;(t), (78) 
"^""^^^ - -^a(t)-mao^(t) + y7^ain(t), (79) 



dt 2 



dt 



= 0, (80) 



^ = 2,{||[Z. + Z(.)]^;(.), (81) 

^ = _2,M[Z., + Z(^)]^;(^). (82) 

In the bad cavity limit where 7c ^ uJm, (r/''^), i^Z [i.e., set {dd{t)/dt) = in 

(79)], the dynamics of the field quadrature, a^{t) + d{t), adiabatically follows 
that of the cantilever position: 

at(i) + a{t) = -i^^Z{t) + -^[di^{t) + aL(i)]. (83) 

7c xTfc 
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Thus monitoring this field quadrature of the cavity mode via homodyne 
measurement corresponds to a measurement of the cantilever position and 
hence the state of the spin. Using the usual input-output relation, [28, 29] 



aout(^) = V^a{t) - Oin, 
we may define an operator corresponding to the output current 
Iout{t) = A/7^/9[aout(t) + aLt(^)] 

t 



(84) 



(85) 



Equation (85) is similar to (42) with ed = 1. By substituting (83) into (85), 
the resultant output current in the bad cavity limit is given by 



(86) 



This equation is also similar to (57) with = 1, obtained from master 
equation approach. 

The Langevin equations for S'^ and S'y effectively decouple from the other 
equations, since they do not appear on the right-hand-side of the equations 
for the other variables. Because of this, they have no effect in our estimate 
of the signal-to-noise ratio, and we shall drop them henceforth. Taking a 
Fourier transform of the linearized Langevin equations, we find, from (85), 
the Fourier component of the output current as 



4ut('^) 



Pic 



iuj -\- 



7c 



m{uj^ — Lip- — iVuj/m) [{iuj — 7c/2) 



[«L(^) - «in(^)] 



(87) 



where G{uj) is the Fourier transform of G{t) = 2rif{t)/\{t). The Fourier 
component of the mean output current signal is then given by 



m 



(88) 
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where 

ic\ ( 2 2 • r 



D{uj) = [lu - f) (ui -u'- ZU-) . (89) 
The output current noise power density spectrum is defined as 

-5out(c^) = y dTe''^^{Ltit)Ltit + T))Git)=o]^ 

= ^ {/ du'e-'^^+^'^\io^t{oo')Iout{oo))GH=o]^ , (90) 

where the subscript G{t) = means evaluation in the absence of the external 
driving force from the spin and {•••}* denotes the time average over t. To 
calculate this noise spectrum, the Fourier transform of the noise correlation 
functions (72)-(75) is needed and given by 

(ain(u;)aL(a;')) = 2nS{u; + u;'), (91) 



1 + coth 



6{uj + Lj'), (92) 



where in obtaining (92) the infinite frequency cutoff limit of the Ohmic ther- 
mal reservoir spectrum, Q ^ oo, has been assumed. After some calculations, 
one can then obtain the output noise spectrum as 




, m J [(7e/2)2+a;2]|/)(^)|2 

^^-%^1^^-thf^H. (93) 
m2 J \D{u})\^ \2kBT J j ^ ' 

The first term in (93), independent of frequency, is the contribution from the 
shot noise of the photons. The next term is the "back-action" noise on the 
position of the cantilever by the radiation (photons). This back action is due 
to the random way in which photons bounce off the cantilever. The final term 
is the thermal noise, due to the thermal Brownian-motion fluctuation of the 
cantilever. Equation (93) is valid at all temperatures. The assumptions made 
in its derivation arc the linearization around the semi-classical steady state 
and the infinite frequency cutoff Vt ^ oo. The high (or medium) temperature 
limit fiuim -C ksT can be obtained by approximating 

Uucoih { ^ 2kBT + (94) 
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We plot these three contributions to the noise for the simulation parameters 
given in section 7. We see that at the oscillator resonance utm = 1, thermal 
noise dominates. 
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Figure 7: We plot the various terms of So^ti^) vs. u, using the parameters 
of section 7. We have scaled 7c/3 = 1. Note that at = uJm = 1 the thermal 
noise dominates. 

Let us define the signal-to-noise ratio per root Hertz as 



We are interested in evaluating SNR(u;) at the frequency equal to the can- 
tilever vibration frequency, u = ujm- Note that 



where the quality factor Q = mUm/T. As a result, the mean output current 
signal (88) at u = Um is enhanced by a factor of (57c/[(7c/2)^ -t- uj'^Y^'^ as 
compared with the u = case. However, a similar enhancement occurs in 
the back-action noise and the thermal noise terms. In other words, driv- 
ing the cantilever at = Um amplifies not only its vibration amplitudes 
due to the the driving force, but also the noise amplitude due to the back- 
action radiation pressure and thermal Brownian motion (see Fig. 7). We find 




(95) 




(96) 
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SNR(cu 



LUrn) can be written as 



SNR(^^) 



(97) 



where 




Hz, corresponding respectively to the spin in the two different states in the 
rotating frame. 

Because the driving force f{t) is periodic, G{u!) is equal to a sum of delta 
functions at cj = cjm, 3cJm, Scj^, . . .. Averaging over a small interval about 
LUra, we Can integrate over the delta function to get a value (for our simulation 
parameters) of SNR(a;^) 220s~^/^. Thus, given a bandwidth of about IHz, 
this should be easily detectable by our measurement scheme. As mentioned 
in section 7, we have assumed a magnetic field gradient roughly two orders 
of magnitude greater than current experiments, and a much lower temper- 
ature. A single spin, therefore, would be below the edge of detectabihty by 
current experimental techniques. Steady improvement in the field strength, 
temperature and spring constant of these experiments, however, should soon 
make single-spin measurement possible. 

If the dominant noise source in MRFM comes from the thermal Brown- 
ian motion of the cantilever, we can estimate the minimum detectable force 
(when the signal-to- noise ratio is one) by keeping only the last term of (98). 
In this case, with a measurement bandwith Au, we obtain from (97), (98) 
and (94) the usual expression of the minimum detectable force at the high- 
temperature limit {hUm <^ ksT) as 



where k = muj'^ is the spring constant of the cantilever. We see, then, that 
improvement can come either from raising the force (by increasing the field 
gradient), lowering the temperature, or lowering the spring constant. 




(99) 
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9 Conclusion 



We have derived an approximate description of single-spin measurement by 
magnetic resonance force microscopy, including both thermal noise and mea- 
surement back-action, and used it to produce numerical simulations of a 
single-spin measurement. These simulations use the quantum trajectory 
method for open quantum systems. The parameters we assumed for this 
simulation were somewhat optimistic; but given the steady improvement in 
experimental technique, we believe that measurements of this type will be 
possible in the near future. 

Single-spin measurements would be very useful in the construction of 
solid-state quantum computers, in which the spin of an electron represents 
a single qubit of information. Given the great interest in solid-state imple- 
mentations as a possibly scalable realization of quantum computers, finding 
practical ways to measure single spins would be very useful. The results of 
our simulations suggest that magnetic resonance force microscopy is a very 
promising approach to this difficult problem. 
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